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Abstract — Using  regular  perturbation  analysis,  we  investigate 
the  propagation  of  a  time-harmonic  acoustic  signal,  generated 
by  a  sinusoidal  boundary  condition,  in  a  half-space  filled  with 
a  classical  thermoviscous  fluid.  It  is  assumed  that  the  flow  is 
described  by  a  recently  introduced,  weakly  nonlinear  partial 
differential  equation  (PDE)  that,  unlike  earlier  models,  exhibits 
a  Hamiltonian  structure  in  the  lossless  limit. 

I.  Introduction 

In  their  recent  article  on  acoustic  propagation  in  classical 
thermoviscous  fluids  (i.e.,  Newtonian  fluids  for  which  the 
heat  flux  vector  obeys  Fourier’s  law),  Rasmussen  et  al.  [1] 
presented  the  latest  addition  to  the  family  of  equations  used 
to  describe  acoustic  phenomena  under  the  finite-amplitude 
approximation.  Assuming  propagation  in  ID  along  the  x- 
axis,  this  weakly  nonlinear  PDE,  which  we  shall  refer  to 
as  the  Rasmussen-Sprensen-Gaididei-Christiansen  (or  RSGC) 
equation,  assumes  the  form 

(c0 ~4>t)4>xx~(t>tt  +  S(t>txx  =  dt[(4>xy  +  co  2(/3—  §)(^t)2]-  (1) 

Here,  cj>  =  <f>(x,  t)  is  the  velocity  potential,  where  u  =  <px 
and  the  velocity  vector  is  of  the  form  v  =  (u(x, f), 0, 0); 
the  positive  constants  Co  and  6  denote  the  sound  speed  in  the 
undisturbed  fluid  and  the  diffusivity  of  sound  [2],  respectively; 
and  /3(>  1),  the  coefficient  of  nonlinearity1,  is  given  by  /3  = 
l  +  f?/(2A),  where  the  ratio  B/A  is  known  as  the  nonlinearity 
parameter  [3], 

Unlike  its  better  known  counterpart  Kuznetsov’s  equa¬ 
tion  [4],  which  in  the  present  context  takes  the  form 

Cq4>xx  ~  </>it  +  Sfttxx  =  dt[(4>x)2  +  Cq  2(/3  —  l)(0t)2],  (2) 

(1)  admits  a  Hamiltonian  structure  in  the  limit  6  — >  0  (i.e.,  the 
lossless  limit),  as  do  the  Euler  equations  to  which  this  limiting 
case  corresponds.  Of  course,  because  so  much  is  known  about 
Hamiltonian  systems,  and  the  fact  that  they  are  amiable  to 
treatment  by  an  array  of  analytical  methods,  a  Hamiltonian 
structure  is  a  highly  desirable  property  for  one’s  mathematical 
model  to  possess.  On  the  other  hand,  Kuznetsov’s  equation  and 
its  variants,  which  have  been  the  subject  of  intense  study  for 
almost  forty  years,  exhibit  a  number  of  interesting  features; 


see,  e.g.,  [5].  As  a  case  in  point,  the  well  known  Burgers’ 
equation  is  readily  derivable  from  (2)  by  assuming,  in  part, 
unidirectional  plane  wave  flow;  see  [4],  [5], 

While  clear  differences  exist  between  (1)  and  (2),  it  is 
important  to  point  out  that  both  are  derived  from  the  usual 
mass,  momentum  (i.e.,  Navier-Stokes),  and  energy  conser¬ 
vation  equations,  augmented  with  the  non-isentropic,  Taylor 
series-based  equation  of  state  [6] 
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B 
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(3) 


for  a  classical  thermoviscous  fluid  undergoing  irrotational, 
compressible  flow.  As  noted  in  [1],  however,  the  lossless 
version  (1)  is  also  derivable  from  a  variational  approach, 
whereas  the  lossless  version  of  (2)  is  not.  Here,  P  is  the 
acoustic  (or  relative)  pressure,  s  is  the  condensation,  and  the 
positive  constant  k  denotes  the  thermal  diffusivity.  It  should 
also  be  noted  that,  in  spite  of  the  fact  that  there  is  no  exact, 
general  equation  of  state  for  liquids,  (3)  is  applicable  to  both 
gases  and  liquids.  The  derivations  of  (1)  and  (2)  are  based  on 
both  the  assumption  that  the  relative  perturbations  about  the 
equilibrium  state  are  small,  but  finite,  and  the  fact  that  (3)  is  an 
analytically  tractable  expression  for  P  that,  while  approximate, 
is  still  able  to  capture  the  essential  nonlinear  physics  involved. 

Another  PDE  that  we  shall  encounter  in  our  investigation 
is  known  as  Stokes’  equation  [7]. 


C0</>: 


O'rxx 


<t>tt  A  ^4*txx  — 


(4) 


which  is  the  linearized  version  of  both  (1)  and  (2).  It  is  of 
interest  to  note  that  this  equation,  which  is  readily  derived 
from  the  linearized  mass,  momentum,  and  energy  conservation 
equations  by  neglecting  the  s2  term  in  (3),  also  describes  the 
transverse  vibrations  of  an  internally  damped  string  [8]  and 
certain  magnetohydrodynamic  (or  MHD)  flows  [9], 

In  the  present  article,  we  carry  out  a  perturbation  analysis 
of  (1)  and  (2)  in  the  half-space  x  >  0  subject  to  the  boundary 
conditions  (BC)s 


<fix(0,t)  =  (70sin(flt),  c/)x(oo,t)=0  {t  >  0).  (5) 


'In  the  case  of  gases,  the  coefficient  of  nonlinearity  can  also  be  expressed  *he  boundary- value  problem  (BVP)  we  consider,  which 
as  /3  =  (7  +  l)/2,  where  y(>  1)  denotes  the  adiabatic  index  [3].  can  be  regarded  as  the  weakly  nonlinear  acoustic  version  of 
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Stokes’  second  problem,  the  positive  constants  Uq  and  O  de¬ 
note  the  amplitude  and  angular  frequency,  respectively,  of  the 
time-harmonic  excitation  applied  at  the  boundary  x  =  0.  Here, 
it  should  be  noted  that  similar  perturbation  studies  involving 
some  of  the  various  other  PDEs  of  nonlinear  acoustics  can 
be  found  in  the  literature;  see,  e.g.,  [10]— [13].  In  particular, 
the  analytical  approach  taken  here  closely  follows  those  of 
Zabolotskaya  et  al.  [12]  and  Jordan  [13],  who  considered  the 
modified  Burgers’  equation  and  the  Darcy-Jordan  poroacous- 
tic  model,  respectively. 

The  primary  aim  of  this  investigation  is  to  compare  the 
RSGC  equation  with  Kuznetsov’s  equation  in  the  context  of 
the  above-mentioned  BVP.  That  is,  we  seek  to  understand 
how  the  differences  in  these  two  PDEs  are  manifested  in 
their  perturbation  solutions.  In  particular,  we  compare/contrast 
the  lowest-order  correction  terms  in  both  the  low-  and  high- 
frequency  regimes.  And,  of  course,  the  exact  solution  of 
(4),  which  is  also  the  zeroth-order  term  in  both  perturbation 
solutions,  is  derived  and  analyzed. 

To  this  end,  the  present  article  is  arranged  as  follows.  In 
Sect.  II,  the  two-term  perturbation  solutions  of  these  equations, 
subject  to  the  above  BCs,  are  determined.  In  Sect.  Ill,  low- 
and  high-frequency  asymptotic  expressions  are  derived  and 
numerical  results  are  presented.  Finally,  in  Sect.  IV,  a  summary 
is  given  and  conclusions  are  stated 

II.  Stokes’  second  problem:  Nonlinear  acoustic 

VERSION 

Re-expressing  (1),  (2),  and  the  BCs  given  in  (5)  in  terms  of 
the  following  nondimensional  variables:  4>  =  cf>/(LUo),  u  = 
u/Uq,  x  =  x/ L,  and  t  =  t(co/L ),  where  the  positive  constant 
L  is  a  characteristic  length,  and  switching  to  a  more  compact 
form  of  notation,  our  BVP  becomes 

□  </>  +  A  (f>txx  =  ^W4t4>xx  +  [{fix)2  +  /(cr)  (<^>t )  2]  t } , 

(x,t)  £  (0,oo)  x  (0,  oo);  (6a) 
4>x(0,t)  =  sin(ujt),  <f)x(oo,t)  =  0,  (t  >  0);  (6b) 

where  n2  :=  dxx  —  dtt  denotes  the  ID  d’Alembertian  operator, 
e  =  Uq/cq  is  the  Mach  number,  to  =  QL/cq,  and  A  denotes 
the  dimensionless  diffusivity  of  sound.  In  addition,  the  index 

fl,  RSGC  equation, 

a:=  {  (7) 

I  0,  Kuznetsov’s  equation, 

and  the  function  f(a)  =  —(cr  —  B / A) /2  have  been  introduced 
here  for  convenience;  0<e<l  follows  from  the  weakly  non¬ 
linear  approximation;  we  observe  that  A  =  1  /Recj,  where  Red 
denotes  the  Reynolds  number  based  on  5;  and  all  superposed 
bars  have  been  suppressed  but  remain  understood. 

With  e  as  the  natural  perturbation  parameter,  and  recalling 
that  both  (1)  and  (2)  were  derived  under  the  finite  amplitude 
approximation,  we  assume  a  regular  expansion  [11],  [14],  [15] 
correct  only  through  O(e);  specifically,  subject  to  e|^i|  <C 
|<)>o|,  we  set  </>  =  where 

:=  4>0  +  e4>i- 


Substituting  (8)  into  BVP  (6)  yields,  after  expanding  and 
equating  like  powers  of  e,  the  following  sequence  of  linear 
PDEs: 


^  4n  T  A  dtxxftn 

°’  ,  n  =  °’  (9) 

which  are  to  be  solved  subject  to 

3*0n(O,f)=  rmM)’  n  ^  dx(j)n{cc,t)  =  0,  (10) 

(0,  n=  1, 

successively  for  n  =  0,1. 

Using  first  the  substitution  </>o(x,  t)  =  ^iXo(x)  exp(— i uit,)  + 
c.c.,  where  “c.c.”  denotes  the  complex  conjugate  of  the  preced¬ 
ing  term,  the  zeroth  order  solution,  which,  as  noted  earlier,  is 
the  exact  solution  of  Stokes’  equation  in  the  setting  of  BVP  (6), 
is  easily  found  to  be 

4>o  =  -|iP0(w)exp[-(a0  -  i/?o)z]  exp(-kuf)  +c.c.,  (11) 

where  Po(oj)  =  («o  —  i/3o)_1  and  the  corresponding  attenua¬ 
tion  coefficient  and  wavenumber  are  respectively  given  by 

-1  +  VI  +  A2w2  _ 

2(1  +  A 2lo2)  ’  ^  ~  * 

(12) 

In  turn,  the  first  perturbation  of  <pQ  is  determined  by  first  setting 

(f>i(x,t)  =  |iXi(x)  exp(— 2iojt)  +  c.c.  and  then  solving  the 
resulting  inhomogeneous  ODE 

Xi(1  —  2\\u>)+Auj~Xi  =  uj[(l  +  ^cr)(ao— ifto)2 f  (ct)]Xq  . 

(13) 

where  a  prime  denotes  A/ Ax.  Omitting  the  details,  it  is  readily 
show  that 

4*1  =  5iPi(w){2exp[-(ai  -  i/3i)x]  -  (ai  -  i/?i) 

xPoM  exp[-2(a0  -  ipo)x)}e~2luJt  +  c.c.,  (14) 


1  +  VUAV 

2(1  +  \2u}2) 


where 


(Aw  +  i)(a0  -i/?0)[l+  2P${u})f(a)] 

PiM  =  - 


and 


ot\  =  2oj\ 


(3i  =  2u\ 


4Aw2(ai  —  i/?i) 

^-1  + V/lT4A2w2 
2(l  +  4A2w2) 


'  1  +  Vl  +  4A2w2 


(15) 

(16) 

(17) 


2(1  +  4A2w2) 

Finally,  using  the  defining  relation  u  =  dx0,  we  find  that 

u^(x,  t)  :=  uo(x,t)  +  eui(x,t),  (18) 


where 


(8) 


u o  =  dx4> o  =  —  exp(— aox)  sin(/3oa;  —  ut)  (19) 


and 


— e 


C.  High-frequency  results 

Once  again  employing  the  binomial  theorem,  but  now 
under  the  assumption  u>  1/A,  the  corresponding  large-w 
0  sin(2/3otc  —  2 wt  +  t)> i)].  (20)  approximations  are  found  to  be 


ui  =  dxfi  =  7?i(w)[e  aiX  sin(/3ia:  —  2u>t  +  ipi) 


Here, 


i?i(w)  =  \Ai  +  \Bi\  and  t/h  =  Arg(Ai  (21) 

where  the  real  and  imaginary  parts  of  Hi+i£?i  =  2Pi(w)(cti  — 
i/?i )  are  given  by 


Ai  = 
and 

B i  = 


(1  +  icr)(a0Aw  +  fa)  -  f(a)(a0Xuj  -  p0)Vl  +  A2w2 

2XJ2 


(1  +  \a)(a0  -  /30Aw)  -  f(a)(a0  +  ftAwj/l  +  AV 

2A^  ’ 


as  ui 

that 


Pn 


A  (n  +  1  )u> 


A  (n  +  l)w 


1  - 

1  + 


1 


2A(n  +  l)u> 
1 


(25) 


2A  (n  +  l)tuj 

oo,  for  n  =  0, 1.  Thus,  from  (22)  and  (25)  we  find 


Vn 


2io 


1  - 


1 


A(n+1)  [  2A(n  +  l)cuJ 

1+  1 


r 

2 

A  (n 

+ 1  w 

2A  (n  +  1  )u 


(26) 


and  where  |  •  |  and  Arg(  • )  denote  the  modulus  and  the  principal  as  to  — >  oo,  for  n  =  0, 1. 
value  of  the  argument,  respectively,  of  a  complex  quantity.  . 
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III.  Analytical  and  numerical  results 
A.  Phase  speeds  and  penetration  depths 

Along  with  the  attenuation  coefficients  and  wave  numbers, 
two  other  quantities  are  of  importance  in  the  study  of  waves 
produced  by  harmonic  excitations.  Specifically,  the  phase 
speeds  and  penetration  depths,  which  are  respectively  given 
here  by 

Vn=  and  dn  =  —  (n=  0, 1).  (22) 

Pn 

These  quantities  are  used  to  characterize  the  structure  and 
behavior  of  the  wave  field  over  space  and  time.  (For  a 
discussion  of  the  physical  significance  of  these  quantities,  see 
[16],  [17],  respectively.) 

In  the  next  two  subsection,  we  derive  low-  and  high- 
frequency  asymptotic  expressions  for  ao,i,  /3o,u  and  the 
quantities  given  in  (22). 


Fig.  1  was  generated  based  on  the  solution  of  (13)  using  the 
software  package  Mathematica  (version  5.2).  The  value  of  /3 
taken  corresponds  to  seawater  at  20°C  and  3.5%  salinity  [3], 
while  those  of  e  and  A  were  chosen  based  primarily  on  the 
need  to  produce  clear,  informative  graphs  over  the  frequency 
range  considered.  In  Fig.  1  we  see  that  the  difference 

AmaxM  :=  max{|Xi||CT=0}  -  max{|A^||CT=i}  (27) 

£>0  x>0 

clearly  increases  as  cu  is  increased.  This  behavior  appears 
to  be  due,  at  least  in  part,  to  the  fact  that  the  coefficient 
If  depends  on  both  uj  and  a.  On  the  other  hand,  we  see 
that  as  l v  is  decreased,  the  two  curves  approach  each  other 
and  start  to  coalesce,  eventually  becoming  indistinguishable. 
Note,  however,  that  decreasing  ui  from  5.0  to  0.1  has  also 
caused  both  quantities  on  the  right-hand  side  of  (27)  to  grow 
so  large  that,  for  the  Mach  number  value  taken,  the  primary 
assumption  upon  which  our  perturbation  solutions  are  based 
is  violated  (see  Fig.  l(c,d)). 


B.  Low-frequency  results 

Suppose  that  u>  < C  1/ A.  Then,  using  the  binomial  theorem, 
it  can  be  established  that  the  attenuation  coefficients  and  wave 
numbers  admit  the  small-tu  approximations 

a„  «  lA(n  +  1  )2w2  [l  -  |(n  +  1)2(16  +  A2)w2] , 

Pn  «  (n+  l)w[l  -  |(n+  1)2(16  -  A2)w2],  (23) 

as  u>  — >  0,  for  n  =  0,1.  Thus,  from  (22)  and  (23)  it  follows 
that 

Vn  «  (n  +  1)  1  [l  +  +  1)2(16  —  A2)w2] , 

dn  ~  A (n  +  1)2uj2  t1  +  +  1)2(16  +  A^w2]  ’  (24) 

as  u)  — >  0,  for  n  =  0, 1. 


IV.  Summary  and  conclusions 

In  this  very  brief  study,  we  have  carried  out  a  regular 
perturbation  analysis  of  the  RSGC  equation  in  the  context 
of  the  acoustic  version  of  Stokes’  second  problem.  We  have 
presented  low-  and  high-frequency  expressions,  examined  the 
linearized  problem,  and  compared  the  RSGC  equation  with 
the  classical  Kuznetsov’s  equation  using  both  analytical  and 
numerical  methods.  Based  on  an  analysis  of  these  findings, 
we  report  the  following: 

1)  The  presence  of  the  thermoviscous  (i.e.,  damping)  term 
Sftxx  iit  (1)  and  (2)  was  found  to  have  a  stabilizing 
effect  in  the  sense  that,  while  the  67(e)  term  in  the 
lossless  (i.e.,  A  — >  0)  case  is  secular  (see  [11,  eqn.  (14)]), 
the  same  term  in  (18)  is  not. 

2)  The  solutions  of  both  (1)  and  (2)  assume  a  diffusive 
character  for  large-cu;  see  Sect.  IIIC.  Consequently, 


e|Xi'(jc)|  (a):  w  =  100.0 


el^i'WI  (b):  o)  =  10.0 


e|Xi'WI  (c):  co  =  5.0 


el^i'WI  (d):  co  =  0.1 


Fig.  1.  e|X(|  vs.  x  for  e  =  0.05,  A  =  0.1.  and  /3  =  3.625.  Solid:  rr  1 . 

Broken:  a  =  0. 


these  PDEs  predict  that  information  contained  in  high- 
frequency  harmonic  signals  is,  practically  speaking, 
rapidly  and  irreversibly  lost,  due  to  the  “smoothing 
property”  of  the  diffusion  equation,  as  such  signals 
propagate  through  classical  thermoviscous  fluids. 

3)  Since  Vo,i  are  strictly  increasing  functions  of  u>,  both  (1) 
and  (2)  predict  that  classical  thermoviscous  fluids  exhibit 
anomalous  dispersion  [16],  with  respect  to  harmonic 
acoustic  signals,  where  it  should  be  noted  that  1  >  Vo  > 
V\  >  0  for  uj  >  0.  (Recall:  Vo,i  are  independent  of  rr.) 

4)  The  only  notable  difference  between  (1)  and  (2)  occurs 


at  higher  frequencies.  Specifically,  Z\max(w)  increases 
with  increasing  oj  (at  least  over  the  frequency  range 
considered  in  Fig.  1).  Otherwise,  the  behaviors  of  (1) 
and  (2)  in  the  setting  of  BVP  (6)  appear  to  be  all  but 
identical.  For  example,  our  numerical  simulations  also 
suggest  that,  for  to  sufficiently  small, 

0<  \Xi\\a=0-  |^IU=1  <1  (Vx  >  0).  (28) 

(This,  of  course,  also  suggests  that  Z\max  is  strictly  non¬ 
negative.) 

5)  Fig.  l(b-d)  clearly  illustrates  that,  below  a  certain  fre¬ 
quency,  whose  value  likely  depends  on  A  and/or  3,  both 
max{|X{||(7=o}  and  max{|2f(| L=i}  start  to  increase 

#>0  (E>0 

(possibly  without  bound)  as  to  is  further  decreased. 
Thus,  within  this  (low-)frequency  range,  there  is  an  in¬ 
dependent  upper  bound  that  e  must  satisfy;  otherwise, 
the  primary  assumption  upon  which  both  perturbation 
solutions  are  based  is  violated  (see  Fig.  1(d)). 
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